Bures Metrics for Certain High-Dimensional Quantum Systems 
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Hiibner's formula for the Bures (statistical distance) metric is applied to both a one-parameter 
and a two-parameter series (n — 2, ... ,7) of sets of 2 n x 2" density matrices. In the doubly- 
parameterized series, the sets are comprised of the n-fold tensor products — corresponding to n 
independent, identical quantum systems — of the 2x2 density matrices with real entries. The 
Gaussian curvatures of the corresponding Bures metrics are found to be constants (— ). In the 
second series of 2™ x 2" density matrices analyzed, the sin(//t/-parameterized sets are formed - 
following an earlier study of Krattenthaler and Slater — by averaging with respect to a certain 
Gibbs distribution, the n-fold tensor products of the 2x2 density matrices with complex entries. 
For n = 100, we are able to compute the Bures distance between two arbitrary (not necessarily 
neighboring) density matrices in this particular series, making use of certain eigenvalue formulas 
of Krattenthaler and Slater, together with the knowledge that the 2™ x 2 n density matrices in the 
series commute. 



PACS Numbers 03.65. Bz, 05.30. Ch, 05.70.-a, 02.40.Ky 

Some five years ago, Hiibner (cf. (2)), in an article entitled, "Explicit computation of the Bures distance for 
density matrices," derived a general formula for the Bures or statistical distance [0 (cLb) for n-dimensional density 
£NJ ' matrices (p). It took the form, 
> 
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where dp is the incremental change in p, and Ai is the i-th eigenvalue corresponding to the eigenvector < i\ of p. 

Dittmann [Q has indicated how formula ([!]) can be reexpressed — using the Cayley-Hamilton identity — in terms 
of certain matrix invariants (thus, obviating any need for the eigenvalues and eigenvectors themselves). Slater Q has 
applied Dittmann's formula for the case n = 3 to a set of /owr-parameter 3x3 density matrices. (More recently still 
||, we have relatedly examined the full eig/zi-dimensional convex set of 3 x 3, that is spin-1 density matrices.) In 
an infinite-dimensional (but three-parameter) setting, Twamley [Q has found the Bures metric for squeezed thermal 
states. Slater Q utilized these results to find the corresponding volume element (which turns out to be simply the 
^ I product of a function of the squeeze factor and a function of the inverse temperature, the phase being irrelevant in 
H ' this regard). The volume elements of Bures metrics are of particular interest in that — if normalizable — they can be 
considered to form prior probability distributions (for purposes of Bayesian inference) over the associated quantum 
systems || . Twamley [Q has also suggested that a physical significance can be attributed to the scalar curvature of 
the Bures metric, in providing a parameterization (coordinate) independent measure of the accuracy of estimation of 
a state, given a reference state. (Scalar curvatures of Riemannian metrics on thermodynamic state spaces have been 
a subject of considerable analysis |lO| , |ll[| .) 

In this communication, we make direct use of the formula of Hiibner (^) to obtain the restriction of the Bures metric 
for two different series of sets of 2™ x 2" density matrices of interest. In general, it requires 2 2n — 1 parameters to 
specify a 2™ x 2™ density matrix. In light of the consequent severe computational demands entailed, we will limit our 
attention here to certain quite special density matrices, requiring but one or two parameters for their specification. 
We bring to the reader's attention, however, Hiibner's remark: "The Bures metric is defined on the whole [emphasis 
his] space of density matrices" p. 224]. 

We examine the associated normalized volume elements and the Gaussian curvatures for these restricted metrics. 
We note that in two dimensions — the framework of our first series of analyses — the Gaussian curvature is simply 
equal to the negative of one-half of the scalar curvature [[l2], p. 184]. In one dimension — the framework of our second 
series of analyses — there is no nontrivial notion of [intrinsic] curvature. 

The starting point for our first series of analyses is the convex set of 2 x 2 density matrices having real entries. Its 
members (n — 1) are expressible as 

1 / 1 + r cos r sin \ , . 
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where (r, 9) are polar coordinates (0<r<l, O<0< 2tt), parameterizing the unit disk. An application of formula 
(|l|) to the density matrices ^ yields the Bures metric, 

dsiPreaU Preai + d/wz) 2 = g r Ar 2 + 2g r gdrd9 + g ge d8 2 , (3) 

where 

9r9 = 0, (5) 

and 

r 2 

9ee = -j- (6) 

We observe that the elements of the metric (Q) - (^|) are independent of the angular parameter 8. Normalizing the 
area (the form, of course, which "volume" takes in two dimensions) element (y/det g = y/g r r9ee) over the unit disk, 
we obtain the prior probability distribution, 

p(r,9) = — ^=L=. (7) 
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r 



We have also applied Hiibner's formula ([!]) to the n— fold tensor product of the density matrix with itself — 
corresponding to n identical, independent two-level real quantum systems — for n = 2, . . . , 7, and obtained, in all 
these six cases, results of the form, 

n 

9rr = 4(1 -r*y (8) 
9re = 0, (9) 

and 

2 

TIT 

9ee = —- (10) 

(The eigenvalues and eigenvectors — to be used in (|l]) — of the n-fold products are directly derivable through basic 
rules from those of p rea i itself.) Since g r g has been found to equal zero, for n = 1, . . . , 7, the polar coordinates (r, 6) 
comprise an orthogonal curvilinear (Lame) coordinate system in these instances, and presumably for all n. (The 
system is not isothermal or conformal [[l3]|l4|], however, in that g rr ^ ggg — not even for a single value of r G [0, 1].) 

We have, following the lead of Twamley , computed the Gaussian curvature (K ) — which equals one- half of 
the negative of the scalar curvature |12|, p . 184] — of the Bures metrics reported above. For our orthogonal (polar) 
coordinate system, this takes the form |l5|, p. 105], 

K = - (—( ^ 9rr ) + -( & 9ee )), (11) 

^y/g rr gee 99 ^g rr gee dr y/g rr gee 

For n — 1, . . . , 7, we have that K = ^, that is, the 2" x 2™ density matrices form spaces of constant positive Gaussian 
curvature. (Braunstein and Milburn |l7| have noted that: "there is an overall V^iV improvement in the precision to 
which we may determine the parameter X as we increase the number of identically prepared systems we can make 
measurements upon. This is familiar to us as the typical improvement upon increasing our sample size; here we see 
that it is a general limit to how well we can determine a parameter from quantum systems.") 

We have also conducted a more limited analysis of the 2x2 density matrices with complex entries. These are 
parameterizable using spherical coordinates (r, 9, <p) in the form, 

1 / l + rcos6> r sin 8 cos <fi — if sin 8 cos <fi \ , . 

Pcompiex = 2 I r sin 6» cos + ir sin 8 cos d> l-rcosfl )' ( ' 
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The associated Bures metric, 



dB(pcom P iex, Pcompiex + dp comp i ex ) 2 = g rr dr 2 + 2g r0 drd6 + geedO 2 + 2g r<t> drd(j) + 2g e<t> d9d<i) + g^dcf 2 , (13) 
for n = 1, has elements (cf. formula (15)]), being independent of the longitudinal coordinate <fi, 

9rr = (14) 

r 2 

9ee = -j, (15) 
r 2 sin 9 2 

500 = 1 , (16) 

and 

g r e = 5r0 = 5e0 = 0. (17) 
The volume element of this metric is normalizable to the prior probability distribution, 

g(r,M)= (18) 

over the Bloch sphere jljj of two-level quantum systems, that is, the unit ball in three-space (0 < r < 1, < 9 < 7r, 
< < 27r). The (Ricci) scalar curvature of this metric is equal to -24. We have been able to compute the Bures 
metric for the two-fold tensor products of p CO mpiex with itself, but the metric elements were given by highly involved 
algebraic expressions, which proved difficult to simplify. However, by setting r, 9, ef> to specific values a number of 
times, we obtained results, in all these cases, fully consistent with the proposition that these elements are, in fact, 



simply twice those given by ((14)-([16|). (The scalar curvature would, then, be equal to -12.) This would adhere to the 
pattern noted above (for n = 2, . . . , 7) with the 2x2 real density matrices. Presumably, there exists a demonstrable 
theorem confirming that this rule holds for all n, for both the real and complex 2x2 density matrices and, possibly 
other types of density matrices, as well. (The space of n x n density matrices, for n > 2, "is not a space of constant 
curvature and not even a locally symmetric space, in contrast to what the case of two-dimensional density matrices 
might suggest" ||.) 

In our other series of analyses, we also apply the formula (0) of Hiibner to a series of 2™ x 2" density matrices 
(n = 2, . . . ,7). Rather than the pair of polar coordinates (r,9), as in the first instance, these density matrices are 
parameterized in terms of a single variable (u or alternatively, (3 — as elaborated below) . The intial (universal quantum 
coding) motivation for studying them was presented in an extended paper of Krattenthaler and Slater |fl8f . These 
2™ x 2" density matrices were obtained by averaging (over the Bloch sphere of two-level quantum systems) the n-fold 
tensor products with themselves — corresponding to n independent, identical systems — of the 2x2 complex density 
matrices (|l2|). (An analogous [unpublished] study has also been conducted, using the 2x2 real density matrices, but 
the proofs of certain propositions have turned out — somewhat surprisingly — to be more problematical, involving 
an intricate triple summation, in that [lower-dimensional] context.) The averaging was performed with respect to a 
one-parameter (it) family of probability distributions, 

T(5/2-u)r 2 sm9 



7 r 3 / 2 r(l-M)(l-r 2 ) t 



In [[L9[ it was argued that this family ( |i9| ) could be given a thermodynamic interpretation by using the changes-of- 
variablc, u = 1 — /3 and r = \/T— e~ E . One, then, arrives at a Gibbs distribution of the form, 

P -PE 

f{E;0) = —n{E), (20) 

where the energy E is taken to be the negative of log(l — r 2 ), the density-of-states or structure function, fl(E), to be 
\fl — e~ E , and the partition function to be 
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(21) 



The natural interpretation of the parameter (3 appears to be that of an effective polarization temperature |T^,g2l (cf- 
©[§)■ 

Explicit formulas were reported in |l8| for the eigenvalues and eigenvectors of the 2™ x 2™ matrices (Cn(/?)) averaged 
with respect to (|l9|). It was found that there are only 1 + \ n/2\ distinct eigenvalues. The eigenvectors constructed in 
p8[ formed bases of the 1 + \n/2\ subspaces, but were not orthogonalized within the subspaces (cf. [^3[ pp. 426-427]). 
The eigenvalues can be expressed as 

X - 1 r(3/2 + /3)r(/3 + g)r(l + /3- g + n) 
n ' ? ~ 2«r(/3)r(l + /3 + n/2)r(3/2 + /3 + n/2)' 9 ~ U ' ' ' ' ' ' W i [ > 

with respective multiplicities, 

_ (n-2d+l) 2 fn+1 



(n+l) 



(23) 



The subspace spanned by the m n<q eigenvectors for the eigenvalue \ rhq corresponds to those explicit spin states |23j : 
sec, 7.5.j] with q spins either "up" or "down" (and the other n—q spins, of course, the reverse). The 2 n -dimensional 
Hilbert space can be decomposed into the direct sum of carrier spaces of irreducible representations of SU(2) x S n - 
The multiplicities (|2^) are the dimensions of the corresponding irreps. The q-th subspace consists of the union of 
™>n,q/{n — 2g + 1) copies of irreducible representations of SU(2), each of dimension (n — 2q + 1) or, alternatively, of 
(n — 2q + l) copies of irreps of S n , each of dimension m„, q / (n — 2q + 1). 

The Bures metrics for the six cases (n — 2, ... ,7) analyzed, take the simple form, gpp&P 2 . We have found that gpp 
equals 

(n = 2) (24) 



4/3(2 + /3)(3 + 2/3) 2 



4/3(3 + /3)(3 + 2/3)2' {n ~ 3) ' (25) 



9(145 + 310/3 + 230/3 2 + 72/3 3 + 8/3 4 ) 
4/3(1 + /3)(3 + /3)(4 + /3)(3 + 2/3) 2 (5 + 2/3) 2 ' 

15(185 + 380/3 + 2 70/3 2 + 80/3 3 + 8/3 4 ) 



(n = 4) (26) 



(n = 5) (27) 



4/3(1 + /3)(4 + /3)(5 + /3)(3 + 2/3) 2 (5 + 2/3) 2 ' 
for n = 6, the ratio of 

45(43 2 60 + 143 6 40/3 + 20 1 740/3 2 + 157 1 70/3 3 + 74 361/3 4 + 2 1 864/3 5 + 38 96/3 6 + 384/3 7 + 16/3 8 ) 

to 

4/3(1 + /3)(2 + /3)(4 + /3)(5 + /3)(6 + /3)(3 + 2/3) 2 (5 + 2/3) 2 (7 + 2/3) 2 . 
and for n — 7, the ratio of 

63(61950 + 200025/3 + 273140/3 2 + 206472/3 3 + 94369/3 4 + 26616/3 5 + 4504/3 6 + 416/3 7 + 16/3 8 ) 

to 

4/3(1 + /3)(2 + /3)(5 + /3)(6 + /3)(7 + /3)(3 + 2/3) 2 (5 + 2/3) 2 (7 + 2/3) 2 . 

For this series of computations, we used MATHEMATICA to obtain fully orthonormal sets of eigenvectors. In doing 
so, for n > 4, computational considerations required us to resort to a somewhat indirect approach, not simply making 
use of the Eigensystem command, but rather the NullSpace command, coupled with our knowledge of the actual 
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eigenvalues. (In all cases, however, it was necessary to, additionally, employ the GramSchmidt command on the 
vectors yielded by the Eigensystem or NullSpace command.) Since our singly-parameterized matrices for a given n 
all commute, we could have relied upon Hubner's formula jl], p. 242], 

d B (p,p + dp) 2 ^ty(dp^ 2 ) 2 , (28) 

"which is simply the Hilbert-Schmidt metric, not on the space of density matrices itself, but on the 'space of roots of 
density matrices' rather" |Q. 

The numerators of the results ((|24|), (|25|)) for n = 2, 3 are simply constants. For n > 3, the roots of the numerators 
also have negative real parts lying between -1 and —n. It is clear from immediate inspection that all the roots of the 
denominators for n = 2, . . . , 7, except for 0, are no greater than -1. There are, of course, singularities of gpp at for 
all these n, and at the other (strictly negative) roots of the denominator. 

In Fig. [l| we plot gpp for n = 2, ... ,7. (There is no nontrivial concept of intrinsic curvature for one-dimensional 
metrics. In this regard, however, it is of interest to note that in his study of squeezed thermal states, Twamley (?j 
finds that the "curvature is independent of the 'unitary' parameters r and 8 and only depends on the 'non-unitary' 
parameter [3.) The curve for n = 7 dominates that for n = 6, which, in turn, dominates that for n = 5, . . . All the 
curves are monotonically decreasing with (3. 




The elements of length (^Jgpp) can be normalized over the range (3 € [0, oo] by dividing ^Jgpp by 7r/6 « .523599 
(n=2), tt/4 w .785398 (n = 3), .987405 (n = 4), 1.1533 (n = 5), 1.29428 (n = 6) and 1.42688 (n = 7). 

From p8| , formula (2.12)], that is, formula (E2) above, we know the eigenvalues of the averaged matrices for arbitrary 
n. Since any two averaged matrices (Cn(/3i), Cn(P2)) for distinct values of (3 are known from |]l8| to share the same set 
of eigenvectors — so, C«(/3i) and Cn(/?2) necessarily commute — the eigenvalues themselves are all that is required 
to compute the (in general, nonlocal) Bures distance between the density matrices. (I thank C. Krattenthaler for 
pointing this out.) This can be deduced from the general formula for the Bures distance 0, 



d B (p 1 ,p 2 ) 2 =2-2tr(p 1 1 /2 p 2 p 1 1 /2 )V 2 . 



(29) 



Employing (E9J) for the case n = 2, by setting pi = C2(/?i),P2 = Qifa), we obtain 



Also, for n = 3, 



dB(C 2 (/3i),C2(/3 2 )) 2 = 2- 



ds(C3(/?l),C3(/3 2 )) 5 



(3lfo 



(3 + 2/3!)(3 + 2/3 2 ) 



(3 + 2/30(3 + 2/32) 



'(2- 


f ft) (2 4 


-ft) 


(3 + 


2ft) (3 4 


-2ft)' 


/(3 


+ ft)(3- 


f ft) 



(30) 



(3 + 2/30(3 + 2/32)' 



(31) 
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In Fig. ||, making a more intensive use of the eigenvalue formula (H|), we plot the Bures distance (ds) between Cwo(Pi) 
and CiooCA)- 




0.5 1 



FIG. 2. The Bures distance (ds) between 2 x 2 density matrices, CiooQ3i) and CiooC/fe), as a function of /3i and fa 
(Of course, the value of ds along the line (3\ — @2 is zero.) The function shown in this figure is computed as 

<M<ioo(ft),Cioo(A)) 



\ 



L»/2J 



2- 2 £ 



(1) \(2) 



(32) 



9 =0 



Given the Bures distance, we would, then, be able to obtain, the Bures metric, through the use of the formula [jlj, p. 
2411, 



1 d^ 

gijdp'dp 3 = ^-^[dB{p,p + tdp) 2 \ t=0 



(33) 



Employing such an approach to finding gpp, would avoid having to compute the eigenvectors of the matrices Cn(P)- 
(We, in fact, found the computation of the eigenvectors of the 256 x 256 density matrices, Cs(/3), to be beyond our 
resources.) 

Let us, making use of (p4|), compare the integrated element of length 



/gppd/3 = tan 



— tan 



(71 = 2) 



(34) 



with the Bures distance itself (|3(]) between C2(/3i) and ^2(^2)- In Fig. Bl we plot the absolute value of this function 
(B4h minus the Bures distance — given by (|3C|). 
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FIG. 3. Excess (A) over the Bures distance 
length (p4j) of the Bures metric. 



1° 

for the case n = 2 — of the absolute value of the integrated element of 



In Fig. [|, we display the analogous result for the case n = 3, making use of ( |3l| ) and the relation, 



/gppdp = tan 1 



tan 



4Wi 



(n = 3). 



(35) 




For n > 3, it appeared it would be necessary to employ numerical integration to evaluate ^fgppAj3 to generate 
similar figures. 

It would be of interest to study the question of whether or not the Bures distance between Cn($0 and Cnifa) can 
be achieved for some particular (geodesic) path through the unrestricted (2 2n — l)-dimensional parameter space of the 
2™ x 2" density matrices (cf. Q). 
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